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We show that parametric context-sensitive L-systems with affine geometry interpretation provide a 
succinct description of some of the most fundamental algorithms of geometric modeling of curves. 
Examples include the Lane-Riesenfeld algorithm for generating B-splines, the de Caste ljau algorithm 
for generating Bezier curves, and their extensions to rational curves. Our results generalize the pre- 
viously reported geometric-modeling applications of L-systems, which were limited to subdivision 
curves. 

1 Introduction 

L-systems were conceived by Aristid Lindenmayer as a mathematical formalism for reasoning about 
growing multicellular organisms lfl4ll . They were originally introduced as an extension of cellular au- 
tomata, allowing for the addition and removal of cells during an automaton's operation, but were soon 
rephrased in terms of rewriting systems lfL5l . This rephrasing led to an elegant definition of L-systems 
using notions and notation borrowed from formal language theory. The relation to formal languages col- 
ored early studies of L-systems, which were often focused on their ability to generate different classes of 
languages lfT0ll29l 

Concurrent with these developments, although initially on a smaller scale, L-systems began to be 
used as a mathematical foundation for the computational modeling of plants. The fundamental observa- 
tion was that a growing filament, i.e., a linear arrangement of cells, can be conveniently viewed as a word 
over an alphabet of cell states. This abstraction was generalized to branching structures and components 
larger than individual cells, leading to astonishingly succinct descriptions of growing plants ll24l . 

In retrospect, the success of L-systems in plant modeling can be attributed to the following fac- 
tors ED : 

1. L-systems specify development in terms of temporally and spatially local declarative rewriting 
rules — context-free and context-sensitive productions — which express developmental processes 
in a concise and intuitive manner. 

2. L-systems describe growing structures in terms of topological (neighborhood) relations between 
structure components, which are automatically maintained when components are added to or re- 
moved from a structure. Consequently, the context for production application is always available. 

3. L-systems refer to model components by their type, state and context, rather than a unique name 
or position in the structure. This simplifies the specification of developmental algorithms. In 
particular, indices are no longer needed. 

These features are not only important to the modeling of biological development, but also make L- 
systems well suited to the description and implementation of some non-biological algorithms [22]. Typ- 
ically, they are characterized by the repetitive application of relatively simple rules to discrete structures 
with a changing number of components. Early studies of such algorithms led to the concise L-system 
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description of linear fractals, including classic space-filling curves |[T9ll25ll3Tll33l . The geometric inter- 
pretations of L-systems used in these studies were based on turtle geometry [1] or chain coding 0. 

More recently, L-systems with an interpretation based on affine geometry ||4j [5j [8] have been 
demonstrated to provide a concise description of the subdivision curves used in geometric modeling EOl 
1261 . Similar to the algorithms for generating fractal curves, subdivision algorithms operate on discrete, 
polygonal lines with a changing number and configuration of components. The rules of change can 
be conveniently described in local terms and formalized as L-systems productions. A comparison be- 
tween L-systems and the traditional specification of subdivision algorithms using indexed points and/or 
matrices reveals the advantages of L-systems ll26l . 

Here we expand the range of geometric modeling applications of L-systems and show that L-systems 
allow for the succinct formulation of some of the most fundamental algorithms for geometric modeling of 
curves. These include the Lane-Riesenfeld algorithm for generating B-splines, the de Casteljau algorithm 
for generating Bezier curves, and their extensions to rational curves. Further examples are discussed in 
detail by Shirmohammadi [30], whose work was a stepping stone for the present paper. Compared 
to OOll . we use a topologically more accurate representation of polygons, rooted in the notion of a 
cell complex [18], with both vertices and edges explicitly represented. This results in a more intuitive 
specifications of the algorithms, reflecting their inherent symmetries. 

2 Preliminaries 

The L-systems employed here are context-sensitive parametric L-systems EHim, extended in two di- 
rections: 

• Parameters are not limited to numbers, but may also be compound data structures, for example 
representing points or vectors in two or three dimensions [11, 26 ] ; 

• Productions may be grouped into subsets (tables ESI ), with a control mechanism deciding which 
subset is applicable in each derivation step |[3l l23H34l . 

The key concepts can be summarized as follows. Parametric L-systems operate on parametric words, 
or strings of modules. Each module is a letter from the L-system alphabet V, which may be associated 
with one or more optional parameters. Beginning with an explicitly defined axiom (initial string), an 
L-system generates a developmental sequence of words using a finite set of productions of the form: 

label : left context < strict predecessor > right context : condition — > successor. (1) 

In each derivation step, productions are applied in parallel to all modules in the predecessor string. If 
several productions apply, the first production in the list is chosen. If no production applies, a module is 
rewritten into itself. For example, the L-system with axiom CO and production p\ to p^, given below 

co: A(1.5)£(2.0,3.0)A(4.5)C(1) 

pi : A{x) :x<2->A(2x+l) 

p 2 : A(x):x>2^B{2x+\) ( 

p 3 : A(w) <B(x,y) >A(z) -> A(w+x)A(y + z) 

generates the word A(4)A(3.5)A(7.5)B(10)C(1) in the first derivation step. In addition to proper L- 
systems, in which the strict predecessor is always a single module, we consider pseudo-L-systems lH9h 
in which the strict predecessor may be a nonempty parametric word over the alphabet V. We also consider 
L-systems operating on circular words, in which the first and last module are considered neighbors. 
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Figure 1: Affine combination of points P(vi) and P(v 2 ). 

Parameters may be t/-tuples of numbers representing positions v in a J-dimensional space (26). An 
affine combination v of m positions vi, v 2 , • • • , v m is given by the expression 

v = aivi + a 2 v 2 H h a m v m , (3) 

where the scalar coefficients a,- add up to 1 : 

ai + a 2 H ha m = l. (4) 

The meaning of the affine combination Q is derived from its transformation to the form 

v = vi+a 2 (v2-vi)H ha„(v m -vi), (5) 

which is a well-defined expression of vector algebra. Specifically, for two positions vi and v 2 we obtain: 

v = a\v\ + a 2 v 2 = vi +a 2 (v 2 - vi) = v 2 + ai(vi - v 2 ), (6) 



which means that point P with position v, noted P(v), divides line P(vi)P(v 2 ) in proportion a 2 : OCi 
(Figure [T]). 

3 B -splines 

A simple example of the application of L-systems to geometric modeling is an L-system specification 
of the Chaikin algorithm 0. Given a (closed) control polygon P(vq)P(vi) . . .P(y m ), this algorithm 
produces a smooth (at the limit) curve, whose shape can be thought of as the result of iteratively cutting 
the corners of the control polygon and its descendants (Figure|2]). This process can be succinctly specified 
by a context-sensitive L-system operating on circular words: 

a: P(v l )P(v 2 )...P(v m ) 

p: J P(v / )< J P(v)>P(v r )^ J P(iv / + |v)P(|v + iv r ) {) 

It is often convenient to explicitly represent not only the vertices, but also the edges of the polygons on 
which the Chaikin algorithm operates. Such a representation is provided by the following modification 
of L-system [7j 

(O: P(vi)EP(v 2 )E ...P(v,„)E 

Pi : P(yi) < E > P(v r ) -> P(|v, + ^v r )£P(Jv/ + |v r ) (8) 
P2-- P(v) -> E 

A point P(y) carries all the information needed to visualize it as, for example, a small circle centered at 
v. In contrast, the visualization of the edge E between points P(v/) and P(v r ) requires an interpretation 
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Figure 2: Illustration of Chaikin's algorithm (adapted from [26]). Beginning with a user- specified control 
polygon (a), the algorithm iteratively removes the polygon corners (red triangles in panels b and c). Each 
old vertex is thus replaced by a pair of new vertices, moved ^ of the distance toward the neighboring old 
vertices (arrows in panel b). An approximately smooth curve obtained after 4 derivation steps is shown 
in panel d. 



rule iPTTl [T2l that is applied at the end of the derivation and gathers the information about the line's 
endpoints: 

h E : P(yi) < E > P{v r ) -> Uyi,v T ) (9) 

Module L now carries all the information needed to draw a line from v/ to v r . We assume that production|9] 
complements all L-systems using explicit edge representation (module E). 
A variant of L-system [8] is given below: 

co: P{v l )EP{v 1 )E...P{v m )E 
p: P(vi)<E>P(v r ) -> EP{\vi + \v r )E 



qi: P{vi)<E>P(v r ) -»• P{\vi + \v r ) 
q 2 : P(y) -»• E 



(10) 



The set of productions is partitioned here into two subsets, labeled p and q. Production p inserts a new 
vertex in the middle of each edge, thus subdividing it into two halves. Productions q\ and 172 replace the 
predecessor polygon with a new polygon that has all vertices placed at the midpoints of the predecessor's 
edges. Derivation proceeds in cycles. Each cycle consists of a single application of production p followed 
by n > applications of productions q\ and ^2 (Figure [3]). For n = 1, such a cycle produces the same 
result as one derivation step in L-system [8] (compare the middle row in Figure [3] with Figure [2]). Lane 
and Riesenfeld have shown that the polygons generated by L-system 10 converge to (uniform) B-spline 
curves of degree n + 1 lTT3l (see also ifTTl 1321 ). B-splines are widely used in geometric modeling due 
to their well understood geometric properties, and the relative ease with which diverse curves can be 
defined by interactively positioning the control points. L-system [10] provides a very concise description 
of this class of curves. 



4 The de Casteljau algorithm 

The de Casteljau algorithm is considered "probably the most important algorithm of all of computer- 
aided geometric design" [6, p. 32]. Given an open control polygon P(v®)P(v2) . . -P(v„) with n > 2 ver- 
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Figure 3: Illustration of the Lane-Riesenfeld algorithm for generating B-splines of arbitrary degree n+ 1 



(L-system 10 1. Beginning with the initial polygon (shown in blue), the algorithm proceeds in cycles 



consisting of a single derivation step using production p (result shown in black) followed by n > steps 
using productions q\ and q 2 (results shown in red). Derivations of a linear (n = 0), quadratic (n = 1) and 
cubic B-spline (n = 2) are shown in consecutive rows. 



tices P located at v°, v°, . . . , v°, the algorithm constructs a polygon P(v\)P(v\) . . .P{v l n _ x ) such that each 
point P(v\) subdivides the line segment P(v{)P(v^ +l ) in proportion t:\-t. The number t is a parameter 
ranging between and 1. This process is iterated n — 1 times, ending with a sing le point P(vJ -1 ). The 
locus of points P(v'l~ l ) obtained by the application of this algorithm for all t G [0, 1] is a curve, called the 
Bezier curve of degree n — 1 defined by the control points (or polygon) P^v^P^v®) . . -P(v®) (Figure 
As with B-splines, the importance of Bezier curves stems from their well understood mathematical prop- 
erties and the ease with which their shapes can be specified or changed by manipulating the control 
polygons. 

The above process of finding a curve point corresponding to the parameter value t can be concisely 
expressed using a simple L-system with axiom CO and productions pi,p2- 



to: P{ Vl )P(v 2 )...P(v n ) 

Pl : P(v)>P(v r ) -»• P{{l-t)v + tv r ) 

Pi- P(v) -> £ 



(11) 



The first production replaces each point in the predecessor string with an affine combinations of this 
point and its neighbor to the right, thus capturing the essence of the de Casteljau algorithm. The second 
production erases the last point in the sequence. The same results are obtained by replacing the right- 
context-sensitive production p\ with its left-context- sensitive counterpart: 



P'l 



P(v/)<P(v) 



P((\-t)v, + tv) 



(12) 



Unfortunately, neither production p\ nor p\ captures the inherent symmetry of the de Casteljau construc- 
tion. As in the case of L-systems[8]and 10 this shortcoming can be addressed by representing the control 
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Figure 4: Illustration of the de Casteljau algorithm (L-system 13 1, a-c) Points obtained for parameter 
values t = 0.2,0.3 and 0.4, respectively. The initial (open) control polygon is shown in blue, the polygons 
obtained in subsequent derivation steps are black, and the resulting point is shown in red. d) The quartic 
Bezier curve obtained as a locus of points generated by the de Casteljau algorithm with the parameter t 
increasing from to 1 in increments At = 0.01. 



polygon as a one-dimensional cell complex: a sequence of points P connected by edges E: 



CO: P(y 1 )EP(y 2 )E...EP(y n ) 

pi: Piyi) < E > P(v r ) -> P((l 

p 2 : E < P(y) >E -> E 

P3 : P(v) -> £ 



t)vi+tv r ) 



(13) 



Production p\ now replaces an edge E with a vertex P dividing this edge in proportions t:\-t. Pro- 
duction p2 performs a dual operation, replacing the vertex between two edges with an edge. Finally, 
production pi, erases the first and the last vertex of the predecessor polygon from the successor polygon. 
The operation of this L-system is illustrated in Figure [4j and the corresponding information flow if given 
by the derivation in Figure [5^. 

An alternative to the sequential evaluation of consecutive points approximating a Bezier curve can 
be described as follows. Points at the inclined boundaries of the derivation tree (Figure [5^) define two 
new control polygons with the same number of vertices as the original control polygon. It is known that 
Bezier curves defined by these new polygons subdivide the original curve into two segments lPT3l (see 
also [6, page 34], for example). Furthermore, the union of the new polygons is closer to the curve than 
the original polygon. A Bezier curve can thus be generated by subdivision, i.e., by constructing two 
control polygons that define the same curve as the original polygon, and iterating this process until a 
desired accuracy of approximation of the curve by the union of its control polygons has been reached. 

The derivation tree for L-system [14] implementing one subdivision level is shown in Figure |5j). The 
consecutively established vertices of the new control polygons propagate from one derivation step to the 
next until the resulting string is complete. The algorithm builds new control polygons sequentially, by 
proceeding from the endpoints of the given control polygon inward. The meeting point — the final result 
of the de Casteljau algorithm — is shared by both new polygons. The L-system is given below: 



co : 


P(v h 2)EP(v 2 ,0)E- 


■P(y n - U 0)EP(y 




Pi ■ 


P(vi,S[) < E > P( Vr 


Sr) -> 


*»((!-*: 


P2 ■ 


E<P(v,s)>E 


:s = -> 


E 


P3 : 


E<P(v,s)>E 


:s^0 -»■ 


IP(v,s)I 


PA ■ 


P(y,s) >E 


:s^0 -»• 


P(v,s)I 


P5 ■ 


E<P(v,s) 


:s^0 -»■ 


IP(v,s) 



(14) 



We now explain the operation of this L-system in detail. To produce the derivation tree in Figure [5J5, 
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Figure 5: (a) Example of the derivation tree generated by L-system 13 which implements the de Castel- 
jau algorithm. A quartic Bezier curve in = 5) has been assumed. The L-system operates on polygons 
defined as sequences of vertices (circles) and edges (horizontal lines). The L-system axiom is given 
in the first row. Solid arrows relate the strict predecessor of each production to its successor. Dashed 
lines indicate the flow of contextual information specifying the location of each vertex. The final result 



is a single vertex shown in red. (b) Example of the derivation tree generated by L-system 14 which 
aproximates Bezier curves of arbitrary degree using repetitive subdivision. The first subdivision level of 
a quartic curve is shown. The central part of the subdivision tree (within the shaded inverted triangle) 
implements the de Casteljau algorithm as shown in panel (a). The left and right parts of the derivation 
tree propagate the vertices and edges established in the previous steps. In order to operate properly, 



L-system 14 distinguishes between vertices and edges of different states and types. For vertices, green 



denotes state s = 0; black: s = 1; red: s = 2. For edges, green denotes type E; black: type /. 

each point is characterized by its position v (first parameter) and state s (second parameter). The states 
represent the following information: 

s=2: an endpoint of a control polygon; 

s=l : a previously established interior vertex of the polygon being constructed; 
s=0: any other point. 

In addition, a distinction is made between edges / that connect pairs of resultant points (including the 
endpoints) and edges E that are still subject to changes by the algorithm. The axiom CO specifies the 
initial control polygon, distinguishing between its endpoints and all other points. Production p\ and 



P2 correspond to productions with the same labels in L-system 13 and capture the essence of the de 
Casteljau algorithm. A new element is function / in production p\, which assigns a state to the resultant 
point. The following rules apply to a vertex in state 0: 

• a vertex adjacent to an endpoint or one interior point becomes an interior point (s=l); 

• a vertex adjacent to two interior points becomes a new endpoint (s = 2). 

The remaining points retain their states. With the assumed coding of states, function / can be written as 

f(st,s r ) = min(st,\)+min(s r ,\). (15) 



The assignment of states allows the L-system to distinguish between points and edges that are still subject 
of the Casteljau algorithm (productions p\ and pi) and points that are ready to be propagated to the 
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Figure 6: Generation of a quartic Bezier curve with the de Casteljau subdivision algorithm, a) The 



initial polygon, b-e) The first subdivision cycle using production set 14 f) Re-initialization preceding 



the second subdivision cycle using production set 16 g-h) The result of the second and third cycle of 
subdivision. Vertex states and edge types are color-coded as in Figure [5] 



resulting string (productions to p$). The propagation itself is effected by the identity productions 
assumed to operate on points and edges to which no other production applies. 

To iterate the subdivision process, the states of points and the types of edges are re-initialized before 
the next subdivision cycle using productions: 



qx : P(v,s) :s = l P(v,0) 
q 2 : I -> E 



(16) 



Given an initial polygon with n + 1 points, the L-system with productions 14 and 16 thus generates a 
Bezier curve of degree n by iterating the cycle of n derivation steps using production set[l4j followed by 



one step using production set 16 The parameter t in production p\ is typically set to 0.5 to achieve a 



relatively uniform distribution of points approximating the curve. An example of a quartic Bezier curve 
generated by this process is shown in Figure [6] 

If the curve degree is known in advance, the result of a complete subdivision cycle applied to a set of 
internal points can be precomputed and encapsulated in a single production. For example, the following 
L-system generates quadratic Bezier curves: 



a: P( Vl )EQ(v 2 )EP(v 3 ) 
p: P(v,)E<Q(v)>EP(v r ) 



0(iv, + \v)EPa Vl + \v + \v r )EQ{\v + ±v r ) 



(17) 



Here we distinguish points by their type (P or Q) rather than state. Points P are the endpoints of a control 
polygon, and are not affected in the subsequent derivation steps. Points Q are the interior vertices and 
can be replaced with another set of vertices in the next derivation step. Generation of a quadratic Bezier 
curve using L-system 17 is illustrated in Figure|7^. 

The principle of this construction carries over to Bezier curves of degrees n > 2. However, the number 
n — 1 of the interior control points that need replacing is then greater than one. Such a replacement can 
be effected most simply using a pseudo-L-system, in which multi-module strict predecessors can be 
rewritten at once |fT9l . For example, a cubic Bezier curve (Figure |7J)) is generated by the following 
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Figure 7: Illustration of the de Casteljau subdivision algorithm, specialized to curves of a given degree. 



a) A quadratic Bezier curve generated using L-system 17 b) A cubic Bezier curve generated using 



L-system[T8l(or 19 1. Points P are shown in red, points Q in black. 



pseudo-L-system: 



CO: P( Vl )EQ(v 2 )EQ(v 3 )EP(v 4 ) 

p : P(yn)E < Q{vi)EQ{v r ) > EP(v rr ) -> 



(18) 



EP \&n + l v l + lv r +lvrr)E 
G(?v, + \v r + \v rr )EQ(±v r + \v rr ) 

An equivalent proper L-system can be obtained by dividing the predecessor and the successor of produc- 
tion p into three parts, for example: 



CO : P( Vl )EQ(v 2 )EQ(v 3 )EP(v 4 ) 

Pi : P(yn)E < Q(vi) > EQ{v r )EP{v r 

P2 : P(v u )EQ( Vl ) <E> Q(v r )EP{v r 

P3 : P(yu)EQ(yi)E < Q{v r ) > EP{v r 



Q{\v n + \vi)EQ{\vu + \vi + \ 
EP{\vu + §Vi + \v r + \v rr )E 
Q{\n + \v r + \v rr )EQ(\v r +\ 



(19) 



Thus, L-systems succinctly express both the basic de Casteljau algorithm and extensions that generate 
Bezier curves of arbitrary or fixed degree using subdivision. 



5 Rational curves 

The algorithms discussed so far have been illustrated with examples of planar curves. However, the 
assumption of planarity is not necessary, and the algorithms operate equally well in three and more 
dimensions. This provides a means of generating curves in space, and also leads to a useful extension of 
Bezier and B-spline curves to their rational counterparts Il6ll27ll. 

Rational curves are generated in a higher-dimensional space, then projected to lower dimensions 
using a perspective projection. Assuming z / 0, the projection of a 3D point P(x,y,z) on the plane z = 1 
from the origin 0(0,0,0) of the underlying coordinate system is the 2D point P'(-, j). Identifying point 
locations with their coordinates, this projection can be accomplished by the interpretation rule 

hp: P(x,y,z) -> PTf/i) (20) 

The subsequently applied edge interpretation rule (production [9]) should operate on the projected points 
P' rather than the original 3D points P, thus taking the form 



h E : P'(v' l )<E>P'(v' r ) -> L(v/,v r ) 



(21) 
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Figure 8: Examples of rational quadratic Bezier curves generated using the de Casteljau algorithm op- 
erating in three dimensions. The endpoints of the control polygon (green) have weight 1. The internal 
control point (red) has weight: a) 0.5, b) 1, and c) 2.5. The weights of the control points are indicated by 
their radius. Example based on [27, pp. 124-125]. 



A 2D point P'(v x ,v y ) can thus be represented by different 3D points of the form P(wv x ,wv y ,w). The 
parameter w is called the weight of point P'. Weights of control points do not affect the shape of the 
control polygon, but bias the resulting curve towards the points with a higher weight. This provides 
an additional means of controlling the shape of rational curves, beyond the manipulation of the control 
polygon. Examples of quadratic rational Bezier curves generated with the de Casteljau algorithm (L- 
system[T3|) using the same control polygon, but different point weights, are presented in Figure[8] 



6 Conclusions 



We have shown that several fundamental algorithms for the geometric modeling of curves can be suc- 
cinctly expressed using L-systems with affine geometry interpretation. Examples include the Lane- 
Riesenfeld algorithm for generating uniform B-splines of arbitrary degree, the de Casteljau algorithm for 
generating Bezier curves, and their extensions to rational curves. Our results generalize the previously 
reported geometric-modeling applications of L-systems, which were limited to subdivision curves [26]. 
Both the previous and current results demonstrate that L-system specifications closely match verbal de- 
scriptions of the modeling algorithms, and capture the information flow underlying their operation. This 
narrows the semantic gap between intuition and mathematical formalism, making L-systems useful as 
a notation for presenting the algorithms. These advantages of L-systems are similar to those observed 
in plant modeling, and stem from the same features. The key feature is the ease of expressing geomet- 
ric algorithms that operate locally on structures with a varying number of components. This ease is 
achieved through index-free notation that emphasizes the topological relations between components (c.f. 
Section [TJ. Whether the use of L-systems could also lead to a notational and conceptual simplification 
of the theorem proofs considered in geometric modeling remains an interesting open question. 
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